Stochastic kinetics of ribosomes: single motor properties and collective behavior 
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Synthesis of protein molecules in a cell are carried out by ribosomes. A ribosome can be regarded 
as a molecular motor which utilizes the input chemical energy to move on a messenger RNA (mRNA) 
track that also serves as a template for the polymerization of the corresponding protein. The forward 
movement, however, is characterized by an alternating sequence of translocation and pause. Using a 
quantitative model, which captures the mechanochemical cycle of an individual ribosome, we derive 
an exact analytical expression for the distribution of its dwell times at the successive positions on 
the mRNA track. Inverse of the average dwell time satisfies a "Michaelis-Menten-like" equation 
and is consistent with the general formula for the average velocity of a molecular motor with an 
unbranched mechano-chemical cycle. Extending this formula appropriately, we also derive the exact 
force-velocity relation for a ribosome. Often many ribosomes simultaneously move on the same 
mRNA track, while each synthesizes a copy of the same protein. We extend the model of a single 
ribosome by incorporating steric exclusion of different individuals on the same track. We draw the 
phase diagram of this model of ribosome traffic in 3-dimensional spaces spanned by experimentally 
controllable parameters. We suggest new experimental tests of our theoretical predictions. 

PACS numbers: 87.16.Ac 89.20.-a 



I. INTRODUCTION 

Ribosome is one of the largest and most complex in- 
tracellular cyclic molecular machines H, [3, H, IJ and it 
plays a crucial role in gene expression [5|. It synthesizes 
a protein molecule, which is a hetero polymer of amino 
acid subunits, using a messenger RNA (mRNA) as the 
corresponding template; this process is called translation 
(of the genetic message). Monomeric subunits of RNA 
are nucleotides and triplets of nucleotides constitute a 
codon. The dictionary of translation relates each type 
of possible codon with one species of amino acid. Thus, 
the sequence of amino acids on a protein is dictated by 
the sequence of codons on the corresponding template 
mRNA. The polymerization of protein takes places in 
three stages which are identified as initiation, elongation 
(of the protein) and termination. In this paper we focus 
almost exclusively on the elongation stage. 

A ribosome is often treated as a molecular motor for 
which the mRNA template also serves as a track. In 
each step it moves forward on its track by one codon 
by consuming chemical fuel [e.g., two guanosine tri- 
phosphate(GTP) molecules]. Simultaneously, in each 
step, it also elongates the protein by adding an amino 
acid; the correct sequence of the amino acids required 
for polymerizing a protein is dictated by the codon se- 
quence on the mRNA template. Therefore, it may be 
more appropriate to regard a ribosome as a mobile work- 
shop that provides a platform for operation of several 
tools in a well coordinated manner. Our main aim is 
to predict the effects of the mechano-chemical cycle of 
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individual ribosomes, in the elongation stage, on their 
experimentally measurable physical properties. We first 
focus on the singie-ribosome properties which character- 
ize their stochastic movement on the track in the absence 
of inter-ribosome interactions. Then we consider the ad- 
ditional effects of the steric interactions of the ribosomes 
and those of the rates of initiation and termination of 
translation on the collective spatio-temporal organization 
of the ribosomes on a track. 

The stochastic forward movement of a ribosome is 
characterized by an alternating sequence of pause and 
translocation. The sum of the durations of a pause and 
the following translocation defines the time of a dwell 
at the corresponding codon. Recently, using an inge- 
nious method, the distribution fit) of the dwell times of 
a ribosome has been measured [6| . We present a system- 
atic derivation of this distribution from a detailed kinetic 
theory of translation which incorporates the mechano- 
chemical cycle of individual ribosomes. 

The exact analytical expression for f{t) which we de- 
rive here is, in general, a superposition of several expo- 
nentials. On the other hand, it has been claimed that 
difference of two exponentials fit the experimentally mea- 
sured f{t) very well. We reconcile these two observations 
by identifying the parameter regime where our theoret- 
ically derived f{t) is, indeed, well approximated by dif- 
ference of two exponentials [7, ^8, 10., JJJ. Moreover, 
we show that {t)~^, inverse of the mean dwell time, satis- 
fies a Michaelis-Menten-like equation p^ . The reason for 
this feature of the mean dwell time is traced to the close 
formal similarity between the mechanochemical cycle of a 
ribosome and the catalytic cycle in the Michaelis-Menten 
theory of enzymes (T2j . 

The elongation of the growing protein by one amino 
acid is coupled to the translocation of the ribosome by 
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one codon. Therefore, {t)~^ is also the average velocity 
{V) of a ribosome on the mRNA track. An analytical 
expression for the average velocity of a molecular mo- 
tor, whose mechano-chemical cycle is unbranched, was 
derived by Fisher and Kolomeisky jlS] in the context 
of motors involved in intracellular transport of cargoes 
. The mechano-chemical cycle of the ribosome in our 
model is, at least formally, a special case of the cycle 
in the Fisher-Kolomeisky model. In this special case, the 
Fisher-Kolomeisky formula for the average velocity of the 
molecular motor, indeed, reduces to the expression for 
(t)^^ in our model of ribosome. 

The average velocity of a ribosome can be reduced also 
by applying an external force (called a load force) that 
opposes the natural movement of the ribosome on its 
track. The force- velocity relation {V){F) (i.e., the varia- 
tion of the average velocity {V) of a motor with increas- 
ing load force F) is one of the most important charac- 
teristics of a molecular motor. Inspired by the recent 
progress in the single-ribosome imaging and manipula- 
tion techniques [1, III [H, [H [H, [13 [13, [HI, we extend 
the formula for {t)~^ appropriately to derive {V){F) for 
single ribosomes. The smallest load force which is just 
adequate to stall a molecular motor on its track is called 
the stall force Fg . We also predict the dependence of Fg 
on the availability of the amino acid monomers and the 
concentration of GTP molecules. 

Our theoretical predictions for f{t), {V){F) and Fg 
show explicitly how these quantities depend on vari- 
ous experimentally controllable parameters. Deep un- 
derstanding of these dependences will also help in con- 
trolling various features of f{t), {V){F) and Fg. In prin- 
ciple, the validity and accuracy of our theoretical pre- 
dictions can be tested by repeating in-vitro experiments 
of ref. [6| for several different concentrations of the amino 
acid monomers and GTP molecules. 

Often many ribosomes simultaneously move on the 
same mRNA track, while each synthesizes separately a 
copy of the same protein. We refer to such collective 
movement of ribosomes on a mRNA strand as ribosome 
traffic because of its superficial similarity with vehicular 
traffic [l^l • In most of the earlier theoretical studies of ri- 
bosome traffic, individual ribosomes have been modelled 
as hard rods and their steric interactions have been cap- 
tured by mutual exclusion [H [H [H, [H, [13, HI [H, Ho, 
[Sll . [3^ . Thus, all those models may be regarded as to- 
tally asymmetric simple exclusion process (TASEP) for 
hard rods [H, [s^] . In some recent works [H, we have 
extended these TASEP-type models of ribosome traffic 
by capturing the essential steps of the mechano-chemical 
cycle of individual ribosomes. We have also reported the 
variation of the average rate of protein synthesis with in- 
creasing population density of the ribosomes on the track. 
In this work we present the phase diagrams of the model 
of ribosome traffic. 

In the earlier TASEP-type models of ribosome traffic 
[H, [H, [H, [M 113, m, [la, [M, [U, m , the phase dia- 
grams were plotted in a two-dimensional plane spanned 



by a and /3, which determine the rates of initiation and 
termination. In this paper we plot the three-dimensional 
phase diagrams of our model of ribosome traffic [s^ in 
spaces spanned by three parameters which, for differ- 
ent diagrams, are selected from a, /3, the availability 
of amino acid monomers and the rate of GTP hydroly- 
sis. Compared to the two-dimensional phase diagram of 
the TASEP-type models of ribosome traffic, these three- 
dimensional phase diagrams provide deeper insight into 
the interplay of single ribosome-mechanochemistry and 
their collective spatio-temporal organization. 

Traffic-like collective movements of ribosomes on a 
mRNA track during translation of a gene was demon- 
strated many years ago by electron microscopy [STj . 
However, to our knowledge, no attempt has been made 
so far to study the phase diagram of ribosome traffic by 
systematic quantitative measurements. But, in contrast 
to most of the earlier works, we have used experimen- 
tally controllable parameters to plot the phase diagrams. 
Therefore, we hope, this paper will stimulate experimen- 
tal studies of the phase diagrams by systematically vary- 
ing the supply of amino acids (monomeric subunits of 
protein) and GTP molecules (fuel of ribosomes) in the 
solution. 

The paper is organized as follows: In section [IT] we 
introduce the model of the mechano-chemical cycle of in- 
dividual ribosomes. The exact dwell time distribution is 
calculated in section Hill while the mean dwell time and 
the physical interpretations of the Michaelis-Menten-like 
equation are presented in section IIVI The connection 
between the mean dwell time and average velocity of a 
ribosome are pointed out in section |Vl where we also 
show the trends of variation of the force- velocity relation 
with variation of some key parameters of the model. The 
variance of the dwell time distribution and the diffusion 
constant of a ribosome are quantitative measures of fluc- 
tuations; the analytical expressions of these quantities 
are presented in section IVII where their relationships are 
pointed out. The distribution of the run times of the ri- 
bosomes on their track and the relation of its first two 
moments with the corresponding moments of the dwell 
time distribution are discussed in section IVIII The ef- 
fects of steric interactions among the ribosomes during 
their traffic-like collective movement on a single mRNA 
track are studied in section IVIIIl the overall rate of pro- 
tein synthesis are presented in subsection IVIII Al while in 
subsection IVIII Bl we plot the 3d phase diagrams of the 
model and also depict 2d projections to compare with 
the corresponding 2d phase diagrams of the TASEP. Fi- 
nally, the results are summarized and main conclusions 
are drawn in section HXl 



II. MODEL OF MECHANO-CHEMICAL CYCLE 
OF RIBOSOME 

The figH] depicts the mechano-chemical cycle of each 
ribosome in the stage of elongation of the protein where 
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FIG. 2: Mechano-chemical cycle of an individual ribosome 
shown in fig[l]is redrawn for the convenience of formulation 
of the master equations. 



constant is fc2 • The peptide bond formation between the 
newly arrived amino acid monomer and the growing pro- 
tein, which leads to the elongation of the protein by one 
amino acid monomer, (and some associated biochemical 
processes, including the arrival of the protein EF-G), is 
captured by the next transition to the state 5 with tran- 
sition rate Ug. All the subsequent processes, including 
the forward translocation of the ribosome by one codon, 
driven by the hydrolysis of another GTP molecule, and 
the exit of a naked tRNA from the ribosome complex are 
captured by a single effective transition from state 5 at 
site j to the state 1 at the site j + 1 with the transi- 
tion rate ujh2- The esssential processes of the cycle are 
summarized in the simplified figure [2l More detailed ex- 
planations of the states and the transitions are given in 
ref.fH. 



III. DWELL TIME DISTRIBUTION FOR A 
SINGLE RIBOSOME: MOST GENERAL CASE 



FIG. 1: (Color online) (a) A cartoon for pictorial depiction of 
the mechano-chemical cycle of an individual ribosome in our 
model. Some of the symbols are explained in (b). 



the integer index j labels the codons on the mRNA track. 
The amino acid monomers are supplied to the ribosome 
in a form in which they form a complex with an adapter 
molecule called tRNA; the complex is called aminoacyl- 
tRNA (aatRNA). Each "charged" aatRNA, bound to 
another protein called EF-Tu, arrives at the ribosome 
from the surrounding medium. The arrival of the cor- 
rect aatRNA-EF-Tu, as dictated by the mRNA template, 
and its recognition by the ribosome located at the site j 
triggers transition from the chemical state 1 to 2 in the 
same location with a transition rate oja- However, if the 
aatRNA does not belong to the correct species, it is re- 
jected, thereby causing the reverse transition from state 
2 to state 1 with transition rate LOp. Hydrolysis of GTP 
drives the transition from state 2 to state 3 with the cor- 
responding rate uJhi- Release of the phosphate group, 
a product of the GTP hydrolysis, is responsible for the 
transition from state 3 to state 4; the corresponding rate 



Because of recent improvements in experimental tech- 
niques, it has become possible to image and manipulate 
single ribosomes d [H [H, 113, E, 13, 113, |2l|. In the 
recent experiments on single- ribosome manipulation 
the distribution of the dwell times of a single ribosome 
at a codon was measured. It was also shown that the 
experimental data fit best to a difference of two expo- 
nentials. More recently, we [111 have demonstrated that 
the numerical data obtained from computer simulations 
of our model can also be fitted to a difference of two ex- 
ponentials. In this section we derive an exact analytical 
formula for the dwell time distribution in our model and 
compare it with the corresponding numerical data ob- 
tained from computer simulations. This analytical for- 
mula shows how the distribution of the dwell times can 
be controlled by tuning the rates of the various sub-steps 
of a mechanochemical cycle of the ribosome. This is a 
new prediction which, in principle, can be tested by re- 
peating the in-vitro single ribosome experiments [6| for 
different concentrations of GTP and aa-tRNA molecules. 

For every ribosome, the dwell time is measured by an 
imaginary "stopwatch" which is reset to zero whenever 
the ribosome reaches the chemical state 1, for the first 
time, after arriving at a new codon (say, j + 1-th codon 
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from the j-th codon). For the convenience of mathe- 
matical formulation, and for later comparison with the 
corresponding results of single molecule enzymatic kinet- 
ics, we make the following assumption: a ribosome finds 
itself in an excited state 1* following the transition from 
the state (j, 5) to (j + 1, 1*) and, then, relaxes to its nor- 
mal state + 1) with a rate constant 6. If the ribosome 
relaxes very rapidly from the state 1* to the state 1, we 
can set 5 ^ oo at the end of the calculation. 

Let P^{j, t) be the probability of finding a ribosome at 
site j, in the chemical state /x at time t. For our calcu- 
lations in this section, we do not need to write the site 
index j explicitly. The time evolution of the probabilities 
Pfiit) are given by 
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dPrM 



dt 



LUgPiit) - UJh2P5{t) 



dPy^jt) 
dt 



i^h2Pb{i) 



(5) 



(6) 



The probability that addition of a new amino acid sub- 
unit to the growing protein is completed between times 
i and t -f Ai is /(t)At. But, 



f(t)M^AP^.{t)^ujn2Pdt)^t 



(7) 



where APi*{t) is the probability that the ribosome is in 
the state 1* in the time interval between t and t + At. 
Therefore, 



dPi.jt) 

= = UJh2Pb{t) 



(8) 



Solving the equations ([T])-®, subject to the normal- 
ization condition 

Pl(i) + P2{t) + Psit) + P4t) + PsW + Pl'it) = 1, (9) 

and the initial conditions 

Pl(0) = 1,P2(0) - P3(0) = P4(0) = P5(0) = Pi.(O) = 0, 

(10) 

we get the time-dependent probabilities P^{t) {fJ. ~ 
1,2,..., 5); the details are given in the appendix. Then, 



FIG. 3: (Color online) Probability density f{t) of the dwell 
times of a single ribosome in the most general case of our 
model for a few different valus of (a) the parameter uja (which 
is proportional to the concentration of tRNA-bound amino 
acid subunits) , and (b) the parameter ujh2 (which determines 
the rate of "stepping"). The continuous curve corresponds to 
the analytically derived expression (lllf) whereas the discrete 
data points have been obtained from computer simulation of 
the same model. 



using the relation ([8]), we obtain the distribution f{t) of 
the dwell times to be 

f{t) = Ciexp{-ujit) + C2exp{-uj2t) 
+ Czexp{—k2t) + Ciexp{-LOgt) + C^exp{—ujh2t) 



where 
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c. = 



UJaU)hlk2(^gUJh2 



{UJI - LOh2){i^2 - LUh2){k2 - UJh2)ii^g ^ i^h2) 



(15) 
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and 



W2 



(w/il 



(17) 



{Uhi +UJp + UJaY 



- ^aOJhl\ 



(18) 

The explicit mathematical formula (|lip for the dwell time 
distribution, which we report in this manuscript, pre- 
dicts how the distribution depends quantitatively on the 
rates of the steps in the mechanochemical cycle of a ribo- 
some. These predictions can be tested by repeating the 
experiments of Wen et al. Q with different concentra- 
tions of amino-acid subunits of the proteins (i.e., aatRNA 
molecules), fuel of ribosome motor (i.e., GTP molecules) 
and ribosomes. 

We plot the distribution pT|) in fi-gO and compare it 
with the corresponding distribution which we have ob- 
tained by direct computer simulation of our model. The 
agreement between the theoretical formula (fTT|) and the 
simulation data is excellent. Note that X]ji=i — 0' 
which implies that f{t) — Q dXt — Q. Moreover, the non- 
monotonic variation of f{t) with t arises from the fact 
that not all of the coefficients are positive. As tUa 
decreases (i.e., effectively, aatRNA become more scarce), 
the tail of the distribution becomes longer and the peak 
shifts to longer dwell times. Moreover, similar trend is 
observed also in the variation of the most probable dwell 
time with the decrease of iOh2- The trend of variation of 
the width of the distribution will be discussed later in 
sectioiiVll of this paper. 



A. Special case 1: ujp = Q 

In the special case ujp = 0, 

f[t) = C^expi-ujat) + C2exp{-uJhit) 
+ C-^exp{-k2t) + C^exp{-ujgt) + C^exp{~ujh2t) 



(19) 



where, is obtained from by replacing loi and 0^2 by 
LHa and w/ii, respectively. The form of the expression (jl9p 
of f{t) makes the underlying physics very transparent- 
f{t) is a superposition of five different terms each of 
which decays exponentially with one of the five rate con- 
stants. Moreover, a clear pattern in the factors in the 
denominators of the coefficients [ji — 1,2,. ..,5) has 
also emerged. 



(a) 



-■+1 COo 

E + S^ li^E* + P 



CO 



E* 



CO 



(b) 



+1 



CO, 



E + S ^ I1 



CO. 



n+1 



-> E* + P 



FIG. 4: (a) Catalytic cycle of an enzyme in the Michaelis- 
Menten theory. E denotes the enzyme while 5* and P denote 
the substrate and product, respectively. The symbol 7i rep- 
resents the intermediate state of molecular complex of which 
the enzyme is a component, (b) Generalization of the cycle 
shown in (a) to n number of intermediate states Ji, /„. 
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FIG. 5: The mechano-chemical cycle of a ribosome, shown in 
the figO is redrawn for the convenience of comparison with 
the MM enzymatic reaction scheme shown in figl^b). The 
symbols /i, 72, /a, 1 4, denote the five intermediate states which 
we labelled in figOby the integers 1, 2, 3, 4, respectively. 



B. Special case II: 



Note that we have derived the general expression (fTTI) 
for f{t) assuming that no two rate constants are equal. 
One can envisage several different possible situations 
where two or more rate constants have identical numer- 
ical values [3^. In order to demonstrate that the form 
of f{t) can get modified under such special conditions, 
in this subsection we consider a very special case where 
u)p = Q and all the nonvanishing rate constants are equal, 
i.e., LOa = ojhi = ^h2 = ^g = k2 = g. lu this case the 
master equations become much simpler and the expres- 
sion for f{t) simplifies to the Gamma distribution 



fit) 



(20) 



where r(fc) is the Gamma function with fc = 5. 
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FIG. 6: The effective mediano-diemical cycie of a ribosome, 
wfiere tlie effective rate constants ijJ2^^ and ujj^/ are given by 
tlie equations ((25} and (|27p . respectiveiy. 



IV. MEAN DWELL TIME: 
MICHAELIS-MENTEN EQUATION? 

Using the expresion for f{t) in 



We get the mean dweh time 
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(21) 
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Further simphfication gives, 
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which is, indeed, the sum of the average time periods 
spent in different steps of the mechano-chemical cycle. 

Next we express the "pseudo" first order rate constant 
uja as LJa = '^ai'^RNA], whcrc [tRNA] is the concentra- 
tion of the tRNA molecules. Then, the eqn. ([23]) can be 
recast as 
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One remarkable feature of the expression is that 
it is very similar to the Michaelis-Menten equation (MM 
equation) for the speed of enzymatic reactions in bulk 
|12j . In chemical kinetics the MM equation is derived for 
the enzymatic cycle shown in figH] where the enzyme E 
enhances the rate of the reaction that converts the sub- 
strates 5* into the product P. In that case the maximum 



rate of the reaction is given by Vmax = while the 
Michaelis constant is (w+2 + 

The steps of the mechano-chemical cycle of an in- 
dividual ribosome, as re-drawn in figEl are very simi- 
lar to those of the generalized MM-like enzymatic cycle 
shown in figHI^b). The fact that the mean dwell time 
for the ribosomes follows a MM-like equation is consis- 
tent with the experimental observations in recent years 
SS SH El El El Ell that the average rate of an enzy- 
matic reaction catalyzed by a single enzyme molecule is, 
most often, given by the same MM equation. 

For our model, we can interpret l/{t) as the average 
rate at which a protein is synthesized by a ribosome, 
where aatRNA plays the role of the substrate and the 
protein elongated by one amino acid is the product. In 
the limit of effectively infinite supply of tRNA molecules, 
on the average, time required to complete one cycle would 
be the sum of the times required to complete the remain- 
ing steps of the cycle each of which has been assumed 
to be completely irreversible. This intuitive expectation 
for the maximum speed of protein synthesis is consis- 



(22) tent with the analytical form 



of 



Further- 



more, in the expression for the Michaelis constant 
the effective rate constants ujj!/ and 1^2''^^ are the coun- 
terparts of oj^i and 0^2 7 respectively, of figlH^a). There- 
fore, so far as the average speed is concerned, the ac- 
tual mechano- cycle, shown in figlSl for a single ribosome 
can be replaced by the simpler MM-like cycle shown in 
figini where uj^ is the counterpart of uj+i, In the limit 
k2 — > 00, — > 00, ujh2 00, the mechano-chemical cycle 
of a ribosome in our model reduces to the enzymatic cycle 
shown in figl3I^a). In this limit, W2^'^ — > ujhi and, hence, 
the expressions (1251) and ((26l) reduce to the correspond- 
ing expressions for Vmax and Km in the MM equation 
for enzymes. 

In reality, however, a ribosome itself is a ribonucleo- 
protein complex that is not an enzyme, but provides a 
platform where several distinct catalysts catalyze the re- 
spective specific reactions. For example, the GTPases 
enhance hydrolysis of GTP molecules while the peptidyl 
transferase catalyzes the formation of the peptide bond 
between the incoming amino acid monomer and the grow- 
ing polypeptide. 



A. Comparison with some earlier works 

Our derivation of the MM-like equation is different 
from the derivation of MM-like equation for cytoskele- 
tal motors reported in ref. [461] where the dwell time 
distribution was not derived. By making one-to-one cor- 
respondence beween the mechano-chemical cycle in their 
generic model for cytoskeletal motors and that in our 
model of ribosome, we find that the MM-like equation 
reported by Bustamante et al. [i^ reduces to the MM- 
like equation 

In a recent work, Jackson et al.fl^l modelled the pro- 
cess of translation as an enzymatic reaction. However, 
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(a) 



FIG. 7: The mechano-chemical cycle of the molecular motor 
in the Fisher-Kolomeisky model for m — 4. The horizontal 
dashed line shows the lattice which represents the track; j and 
j + 1 represent two successive binding sites of the motor. The 
circles labelled by integers denote different "chemical" states 
in between j and j + 1. 
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there are crucial differences between their formulation 
of translation and our interpretation of the mechano- 
chemical cycle in our model. In their formulation, Jack- 
son et al.[47[ treated the completely synthesized protein 
as the product of the enzymatic reaction, i.e., the run of 
a single ribosome from the initiation site to the termina- 
tion site was treated as a single enzymatic reaction. In 
contrast, translocation of a ribosome from one codon to 
the next, and the associated elongation of the growing 
polypeptide by one amino acid has been treated in our 
calculation here as a single enzymatic reaction. 



V. FORCE- VELOCITY RELATION 

Utilizing an earlier result of Derrida |48i] . Fisher and 
Kolomeisky proposed a general formula for the average 
velocity (V) of a generic model of molecular motor where 
the mechano-chemical transitions form unbranched cy- 
cles. Each cycle consists of m intermediate "chemical" 
states in between the successive positions on the track of 
the motor (Figl7|). The forward transitions take place at 
rates Uj whereas the backward transitions occur with the 
rates Wj . Choosing the unit of length to be the separa- 
tion between the successive equispaced positions of the 
motor on the track, the average velocity (V) of the motor 
is given by [ll] 



V = 



1 

Rrr. 



3=0 



(28) 



where 



J=0 



3=0 



rn — 1 k 

En 

fc=l i=l 



(29) 



Formally, our model of ribosome is a special case of the 
Fisher-Kolomeisky model where uq = LOa, ui = u)hi, 1*2 = 
^2, U3 = ujg, tt4 = ujh2 and wi = ujp. Hence, in this 
special case the eon. ([28)1 can be written in a compact 
form as 



V 



(30) 



(b) 
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FIG. 8: (Color online) Force-velocity relation for a ribosome 
in our model for a few different values of the parameter (a) cua 
(which is proportional to the concentration of tRNA-bound 
amino-acid subunits), and (b) uig (which can be controlled by 
varying GTP concentration). The continuous curve has been 
obtained from the formula (|34p whereas the discrete symbols 
denote the numerical data points obtained from computer 
simulations of the model. 



with 



and 



^h2 — <^h2/keff 



UJp 



1 



1 



(31) 



(32) 



Note that fc~^^ is an effective time delay induced by the 
intermediate biochemical steps in between two successive 
hoppings of the ribosome from one codon to the next 
'36] . Interestingly, simplification of the exact expression 
([23)1 yields the same formula ((30|) which we derived as a 
special case of the Fisher-Kolomeisky formula for average 
velocity. 

In our model the load force will only affect the 
mechanochemical transition from state 5 at j to state 
1 at j + 1. The dependence of the rate constant ujh2 on 
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F is given by 



uJh2{F) = uJh2{Q)exp[ - 



F S 



(33) 



where u;^2(0) is the magnitude of the rate constant ujfi2 
in the absence of load force and the typical length of each 
codon is (5 = 3 X 0.34nni. Thus, when subjected to a load 
force F, the force velocity relation for a single ribosome 
becomes 



ViF) = 



^h2{F) 

l + ^h2{F) 



(34) 



The force- velocity relation {V){F) has been plotted in 
figlHI^a) and (b) for a few different values of uja and Wg, 
respectively, to demonstrate the dependence of {V){F) 
on the supply of amino acid monomers and the chemical 
fuel GTP. For fixed uja and Wg, {V) decreases with in- 
creasing F and vanishes a,t F = F^ which is identified as 
the corresponding stall force. Moreover, for a given F, 
ly) increases monotonically with increasing Wa and uig 
although the rate of increase gradually slows down. It 
is interesting to note that Fs is independent of both uja 
and ujg because, at stall, a ribosome uses neither amino 
acid monomers nor GTP. For the typical values of the 
rate constants, which we have used in figlHl ~ 25 — 27 
pN. This theoretical estimate is consistent with the value 
26.5 pN reported by Sinha et al. 
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VI. FLUCTUATIONS: MEAN SQUARE DWELL 
TIME AND DIFFUSION CONSTANT 

A. Fluctuations in dwell times 



FIG. 9: (Color onlne) The randomness parameter r, defined 
by the equation ((39}, is plotted against the concentration of 
aa-tRNA for a few different values of ooh (in (a)) and uip (in 
(b)). 



Mean-square dwell time is defined by 



{t') = / t'f{t)dt. 
Jo 



For our model 

= 2 

LU.2 

The expression can be expressed also as 



C\ C2 ^3 ^4 



where, 



P 



1 



1 



1 



LOa LOhlJ \K2 ^3 ^^h2, 
1111 



UJaUJhl ^a^cj '^a^h2 

1 1 1 



07^1^2 ^hll^g ^hl^h2 
1 1 1 



fczWg k2UJh2 UJgUJh2 



(35) 



(36) 
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(38) 
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FIG. 10: (Color online) The randomness parameter r, defined 
by the equation ((39}, is plotted against the concentration of 
aa-tRNA for a set of parameter values where all the magni- 
tudes of all the rate constants, other than uja, are quite high. 
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Note that only the first term involves ujp. The remain- 
ing ten terms are inverse of the products of the five rate 
constants. 

Let us define "randomness parameter" r as 

(39) 

Note that r is a quantitative measure of the fluctuations 
in the dwell times of a ribosome. By substiuting the 
expressions of (^t'^'j and (t) into we obtain, 



{t') {ty 



{ty 



(40) 



A non-trivial feature of the expression (I40p is that it can- 
not be obtained simply by substituting wj^/ and lo2^^ 
into the expression for r derived by Kou et al. [i^l for the 
two-step Michaelis-Menten enzymatic reaction. In other 
words, the fluctuations of the dwell times in the five-step 
model for the kinetics of ribosomes cannot be captured 
by the effective two-state model drawn in figlH 

The randomness parameter r is plotted in figlQ] as a 
function of the tRNA concentration for a few different 
values of the parameter ojh (in (a)) and ojp (in (b)). We 
find (not shown in any figure) that the numerator of r 
(i.e., {t^) — {ty) decreases monotonically with increasing 
concentration of tRNA; it is the variation of the denom- 
inator of r with tRNA concentration that is responsible 
for the non-monotonic variation of r. 

It is well known [42] that, for a one-step Poisson pro- 
cess, r = 1. At extremely low concentrations of aa-tRNA, 
the binding of a correct species of aa-tRNA to the A site 
on the large subunit of a ribosome is the rate-limiting 
step in its mechano-chemical cycle. Therefore, r is unity 
at sufficiently low values of aa-tRNA. r decreases with 
the increase of aa-tRNA concentration. This decrease 
is caused by the formation of intermediate complexes 
which also affect the rates of progress of the mechano- 
chemical cycle. However, with the further increase of 
aa-tRNA concentration, the randomness parameter r in- 
creases again. Finally, the randomness parameter satu- 
rates to a value which is determined by the number of 
rate-limiting steps in the mechano-chemical cycle. Such 
nonmonotonic variation of r with aa-tRNA concentration 
reduces to a monotonic decrease when the magnitudes of 
the rate constants are sufficiently high (see fig fTO]) . 



B. Diffusion constant 

The diffusion constant is a measure of fiuctuations 
around the directed movement of the ribosome, on the 
average, in space. We now derive a closed form expres- 
sion for D and relate it to the fluctuations in the dwell 
times. Fisher and Kolomeisky's general result for diffu- 
sion coefficient D is 



D 



{VSn + dUn) [N + 2)V 



7?2 



d 

N 



(41) 



where 



and 



while, 



N-l N~l 
j=a k=Q 



N-1 



/ N-1 k 



k=l i=l 



J- J— I 



N-1 
Rn Tj 



(42) 



(43) 



(44) 



(45) 



In our units d — 1. Therefore, in our model the expres- 
sion for D becomes, 

D^{{ty-2^,)/2{ty (46) 

Finally, we observe that r, which is a measure of the 
fiuctuations in the dwell times, is related to D and {V) 
by [5^. 



2D 



(47) 



VII. DISTRIBUTION OF RUN TIMES 

In this section we report the distribution of the run 
times r of an individual from the start codon to the stop 
codon. The run time is related to the dwell times by the 
relation 



L 



(48) 



Central limit theorem states that, as L — > oo, the distri- 
bution G{t) of the run times r approaches a Gaussian, 
irrespective of the nature of the distribution of the dwell 
times, since the dwell times at different codons are inde- 
pendent of each other. Obviously, for sufficiently large L 



G{r) 



1 



v/(27ra2) 



where 



and 



{T)=L{t) 



(49) 



(50) 



(r^) - {ry = L{{t^) {tyy (51) 

Using {t) from and (i^) from ^7^, we obtain the 
Gaussian distribution G{t). The gaussian distribution 
G(t) thus obtained is plotted in figHH it is in excel- 
lent agreement with the corresponding numerical data 
obtained from computer simulations. 
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J, there is another ribosome at site k. Then, Q{i\k) = 
1 — P{i.\k) is the conditional probabiUty that, given a 
ribosome in site j, site k is empty. In the mean- field 
approximation, the Master equations for the probabilities 
P^(j, t) are given by ^36,] 



dPi{j,t) 
dt 



- uJh2Pr.{j-l,t)Q{]A_\3-l+i) 

+ L0j,P2{j,t)-UJaPl{i,t) 



(52) 



dP2{j,t) 
dt 



= LOaPlU, t) - [iOp + UJhl)P2U, t) (53) 



FIG. 11: (Color online) Distribution of the run times of a 
ribosome in our model. The continuous curve is the Gaussian 
distribution predicted by our theory while the discrete data 
points have been obtained from computer simulations. The 
inset shows the dwell time distribution for the same set of 
parameter values. 



VIII. EFFECTS OF STERIC INTERACTIONS 
OF RIBOSOMES 



The average velocity of a ribosome is also the mean 
rate of polymerization of a protein. We define the flux 
of ribosomes to be the total number of ribosomes leaving 
the stop codon (i.e., j — L) per unit time. Obviously, the 
overall rate of protein synthesized from a single mRNA 
template is identical to the flux of the ribosomes on that 
mRNA track. The number density of the ribosomes is 
given by p — N/L. The size of a typical ribosome is such 
that, simultaneously, it covers £ codons where £ ^ 1. 
We treat £ as a parameter of the model. For a given 
number N of ribosomes, the total fraction of the lattice 
covered by all the ribosomes is given by the coverage 
density pcov = N£/L. 

In the preceeding sections, we have ignored the possi- 
bility of steric interactions among the ribosomes. Con- 
sequently, the average velocity was independent of the 
ribosome population on the given mRNA track. Such a 
scenario holds at most at sufficiently low coverage den- 
sities. However, in the presence of inter-ribosome inter- 
actions the average velocity becomes a function of the 
coverage density thereby giving rise to non-trivial varia- 
tion of the flux J (and, hence, the overall rate of protein 
synthesis) with Pcov Moreover, the density profile of the 
ribosomes on the track also exhibits interesting features. 
In this section we study the spatio-temporal organization 
of the ribosomes in terms of the flux as well as the den- 
sity profiles on a single mRNA track and plot the phase 
diagrams of the model. 

Let Pfi{j, t) be the probability of finding a ribosome at 
site j, in the chemical state p at time t. Also, P{j,t) — 
'^^^i Pfi{j, t), is the probability of finding a ribosome at 
site j, irrespective of its chemical state. Let P{i\k) be 
the conditional probability that, given a ribosome at site 



dP3{j,t) 
dt 

dPiU,t) 
dt 



= ^hlP2{],t)~k2P5{3,t) 
= k2P3{j,t)~UJgP4{j,t) 



dP^{],t) 



L0gPiij,t)^L0h2P5U,t)Q{i\j+£) 



dt ' 

Because of the normalization condition 

p{j,t) = j2p,u^t) = ^ = p 

not all of the five P^(j, t) are independent. 



(54) 
(55) 
(56) 

(57) 



A. Effects of steric interactions on rate of protein 
synthesis 

The dwell time distribution f{t) certainly gets affected 
by the steric interactions. As a first step, we have calcu- 
lated the effects of the interactions on the average veloc- 
ity which is just the inverse of the mean dwell time. 

Under periodic boundary conditions, in the steady 
state, Pp(j, t) become independent of j and t. From the 
steady-state limit of the equations ([52ll - (|56|) . we derived 
the expressions P^ {p = 1, 2, ...5) and the flux ^36.] 



JpBC = p{V) 



i^h2p{l ~ pl) 



{I + p - pi) + VLh2{l ~ pi) 



where 



n 



h2 



t^h2/keff. 



(58) 



(59) 



and kef / is given by (l32|) . 

From ([55)1 the p-dependent average velocity (V) can be 
obtained by dividing J by p. At sufRciently low densities 
this expression reduces to 



J = 



ii + n 



h2) 



(60) 



and, hence, we recover the exact formula (|5D|) for the 
average velocity of a single ribosome. 
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I: HD-LD phase boundary 
II: HD-MC phase boundary 
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FIG. 12; (Color online) A Sd-phase diagram of our ribosome 
traffic model. The LD and HD phases coexist on the surface 
I (colour red). Surfaces II (colour purple) and III (colour 
blue) seperate the MC phase from the HD and LD phases, 
respectively. Phase diagram is shown in (a) and (b) from 
two different orientations. The parameters used are Up = 
0.0028s-\a;a = 25.0, = 2As-^ = 25.05"^ 



B. Phase diagrams under open boundary 
conditions 

Initiation and termination of protein synthesis are cap- 
tured more realistically by imposing open boundary con- 
ditions (OBC) than by the periodic boundary conditions 
(PBC). Whenever the first £ sites on the mRNA are va- 
cant this cluster of sites is allowed to be covered by a fresh 
ribosome with the probability a in the time interval At 
(in all our numerical calculations we take At — 0.001 
s). Since a is the probability of initiation in time At, the 
corresponding rate constant (i.e., probability of initiation 
per unit time) uja is related to a by a = 1 — g-'^t^x^*. 
Similarly, whenever the rightmost i sites of the mRNA 
lattice are covered by a ribosome, i.e., the ribosome is 
bound to the L-th codon, the ribosome gets detached 



FIG. 13; (Color online) Same as in fig |12l except that luh has 
been used instead of Ua- 



from the mRNA with probability [3 in the time interval 
At; the corresponding rate constant being denoted by 
Up. For all further discussions in this paper, we'll as- 
sume w/ii = U!h2 = i^h because both these processes are 
driven by GTP hydrolysis. 

TASEP is known to exhibit three dynamical phases, 
namely, high-density (HD) phase, low-density (LD) 
phase and the maximal current (MC) phase in the a — /3 
plane. Our main interest is to explore the nature of the 
dynamical phases in different regions of the 4-dimensional 
space spanned by a, /?, Pui^ , P^j^ ■ 

The parameters Ua and Uh can be controlled by varying 
the concentrations of the aa-tRNA molecules and GTP 
molecules in the solution. The parameter a is determined 
by the rate of assembling of the large and small subunits 
of a ribosome, their final coupling on the initiation site 
and the assistance of several other regulatory proteins in 
the initiation of the actual polymerization of a protein. 
Strictly speaking, a single parameter P captures essen- 
tially two different events both of which take place at 
the termination site j = L. After the full protein has 
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with particle densities p- and p+, respectively. These 
two reservoirs are essentially two infinite lattices with 
the number densities p- and p+, respectively. We calcu- 
late the unknown densities p- and p+, in terms of the 
rate constants of our model, by imposing the requirement 
that these reservoirs give rise to the same probabilities 
a and /3, of hopping with which a ribosome enters and 
exits, respectively, the open system. 

The extremum principle then relates the flux J in the 
open system to the flux J(p) for the corresponding closed 
system (i.e., the system with periodic boundary condi- 
tions). Extremum current hypothesis [H, [H, HI, [s^l 
states that, for the open system connected to the two 
reservoirs of number densities p_ and p+ at its entrance 
and exit, the flux J is related to the corresponding flux 
JpBC in the closed system by 



J 



max JpBcip) P- > P > P+ 
min JpBcip) if P- < p < p+. 



(61) 



Since the flux-density relation (also called the fundamen- 
tal diagram) of our model of ribosome traffic under pe- 
riodic boundary conditions exhibits a single maximum, 
the extremum principle reduces to a simpler maximum 
current principle (MCP). According to this MCP, in the 
limit L — > oo, 

J" = max JpBcip) ii p- > p > p+- (62) 



1. Calculation of p- and p+ 



FIG. 14: (Color online) 2d-cross sections of tlie 3d pliase dia- 
grams of the ribosome traffic model parallel to the a-P plane 
for different values of (a) Pu)^ , and (b) Puj^ , all projected 
onto the same figure. The numbers on the phase boundary 
lines represent different values of (a) Paj„ and (b) Pui^ . The 
HD and LD phases coexist on the curved lines whereas the 
straight lines separate the MC phase from the HD and LD 
phases, as shown. The parameters used for this figure are 
LOr, = 0m28s-\uJa = 25.0, fc2 = 2As-\uJa = 25.03"^ 



been polymerized, the ribosome releases the protein into 
the surrounding medium and then dissociates from the 
mRNA track (the decoupling of the two subunits also 
takes place; these are then recycled for another round of 
protein synthesis) [H, Issl. T herefore, the value (3—1, 
which we assumed in ref . [36| is, in general, not very re- 
alistic. Even in the special case /3 = 1, in ref.[3a| we 
reported only a couple of two-dimensional cross sections 
of the full phase diagram of this model. In this paper we 
plot phase diagrams in three-dimensional spaces spanned 
by a — /3 — P^^^ and a — P — P^^^ . 

For plotting the phase diagram, we use the same ex- 
tremum principle [54, 55, 56, 57] which we used in ref. [3^ . 
In this approach, we imagine that the left and right 
boundaries of the system are connected to two reservoirs 



From (|58p , the maximum flux under PBC corresponds 
to the number density 




1 



V^(i+TM + i 



(63) 



Next we calculate p- . We use symbol 1 to represent the 
sites covered by ribosome while the symbol represents 
the sites which are not covered by any ribosome. Let 
pjump ^j^g probability that, given an empty site, from 
left a ribosome will hop onto it in the next time step. We 
have 



Pi 



Pj l 1 10) X Ps X cj;,2 X A t (64) 



where P5 is the probability of finding ribosome in state 
5 inside the reservoir and the conditional probability 
P{ 1 1 10) represents that, given an empty site, 

e 

there will be a ribosome in the adjacent £ sites to its left. 
On the basis of arguments similar to those presented in 
ref.[36j, we get 



(65) 



(1 + P - P^) 
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Moreover, as argued in ref.[36[, 



1 



Now, p_ is the solution of the equation a = p^"^P and, 
hence, we get, 

where Puj,^^ the probabiUty of hydrolysis in time A t. 

Following similar arguments, we now calculate p^. The 
probability P|"™p^ given that a ribosome which covers £ 
successive sites, will hop onto the next adjacent empty 
site to its right in the next time step, 



Pl""'P = P( l 1 1 0) X Ps X P^,, 



(68) 



where P( l 1 1 0) is the conditional probability of 

finding a hole at site j, given that the site {j — £ — 1) 
is occupied by the leftmost part of the ribosome. It is 
straightforward to show that 



^|0) 



1- pi 



(69) 



i + p-pe 

Now, p+ is the solution of the equation f3 = p^"'-P-^ and 
hence, we get 

P+ 



I3{i + n 



h2) 



P. 



I3{i + nh2){e-i)-ep^,, ^'^^^ 

In the hmit keff — » oo, the expressions (170)) and ([S7|) 
for and p_ reduce to the corresponding expressions 
P- — a and p+ = 1 — /3, respectively, for TASEP. 



2. Surface separating LD and MC phases 



The MCP imposes the condition 

P- = P* 



(71) 



on the surface which separates the LD and MC phases. 
Substituting the expressions (j67p for p_ into equation 
([7T|) we obtain 



P^h2P* 



(72) 



(i + r!,,2)(i-p*(^-i)) 

3. Surface separating HD and MC phases 
From the MCP 

P* = P+ (73) 

on the boundary between the HD and MC phases. Using 
the expression ([70]) for p+, we obtain 



(74) 



4- Surface of co-existence of HD and LD phases 

Since the HD and LD phases co-exist on the surface 
separating these two phases, we obtain the boundary by 
solving the equation 



JpBc{P-) = JpBcip^ 



(75) 



because the same current passes through the two coex- 
isting phases in the steady state, where the density on 
the entry side is p- and that on the exit side is p+ . 

Now incorporating the expressions of p_ from eqn. 
([67]) and p+ from eqn. (f70|) into eqn. ([58]) and using eqn. 
(175)1 we find that the equation of the surface of coexis- 
tence of LD and HD phases is given by 



a 



p^,,Pii + r!,2) 



p^,j + p{i-e + 2nh2 - mh2 + ^l^) 

or, equivalently. 



(76) 



/3 



a £ P^ 



(1 + nh2){Pu^^ -a + £a- VLh2a) 



(77) 



5. Phase diagrams 

In figd^ we have plotted a 3d phase diagram of the 
ribosome traffic model, in the a — (3 — P^j^ space, which 
we obtained by following the MCH-based approach ex- 
plained above. The corresponding 3d phase diagram in 
the a — fi — Pi^f^^ space in plotted in figUHl The LD and 
HD phases coexist on the surface I. A first order phase 
transition takes place across this surface. Surfaces II and 
HI seperate the MC phase from the HD and LD phases, 
respectively. The 3d phase diagrams plotted in figsfTWa) 
and ll3f a) are differently oriented in figs fTW b) andfTWb). 
respectively, to show the regions hidden in fig fTW a) and 
[TST a) behind the surfaces I, II and HI. 

By drawing flat surfaces parallel to the a- (3 plane, each 
corresponding to a fixed value of P^^^ (in (a)) or Pi^^^ (i'^ 
(b)), we have obtained the curves of intersection of this 
flat plane with the surfaces I, II and HI. By projecting 
these curves on the plane Pu,,^^ = Oj ^Iso obtained 
the 2d phase diagram of the system in the a-j3 plane for 
several different values of P^j^^ ■ This phase diagram helps 
in comparing and contrasting our results for the ribosome 
trafhc model with the 2d phase diagram of the TASEP 
in the Pj^^^ plane fFig fH)) . The most interesting feature 
is that, unlike TASEP, the lines on which HD and LD 
phases coexist are curved. This characteristic seems to be 
the general feature of such phases diagrams, rather than 
an exception; similar curved lines of coexistence between 
HD and LD p hases have been observed also in some other 
contexts |58| . 
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f3 > (i+Oh2)(P^^2-a+to-nfc2a) " < (1+J1h2)(l-P.(^-1)) 



ff. Average density profiles 

The bulk density of the system is governed by the fol- 
lowing equations: 




IX. SUMMARY AND CONCLUSION 

In this paper we have derived the exact analytical ex- 
pression for the distribution of the dwell times of ri- 
bosomes at each codon on the mRNA track. For this 
purpose we have used a model that captures the essen- 
tial steps in the mechano-chemical cycle of a ribosome. 
As more details of this cycle get unveiled by new ex- 
periments, our model can be extended to capture those 
new features and the dwell time distribution can be re- 
calculated accordingly. Moreover, some of the transitions 
in the mechano-chemical cycle used in our model may 
require reinterpretation to reconcile with new observa- 
tions. Neverthless, at this stage, the dwell time distri- 
bution predicted by our theory agrees qualitatively with 
the corresponding distribution observed in-vitro single ri- 
bosome experiments. Moreover, our prediction can be 
tested quantitatively by repeating the single ribosome 
experiments varying the supply of amino acid monomers 
and GTP molecules. 

From the full distribution, we have also calculated the 
mean dwell time which satisfies a Michaelis-Menten-like 
equation. We have pointed out the formal similarities 
between the cycles, and the corresponding equations, for 
a single enzyme molecule and a single ribosome, which 
are responsible for the Michalis-Menten-like form of the 
mean-dwell time. The inverse of the mean-dwell time is 
also the average velocity of the ribosome. The expression 
of this average velocity obtained from the dwell time dis- 
tribution is identical to that obtained by an alternative 
approach pioneered by Fisher and Kolomeisky in the con- 
text of generic models of molecular motors. Finally, fol- 
lowing standard procedure, we capture the effects of load 
force by modifying the rate constant U!h2 and predict the 
force- velocity relation and its dependence on experimen- 
tally controllable parameters. From this relation we have 
estimated the stall force of a ribosome. Our theoretical 
estimate is consistent with the experimentally measured 
value reported in the literature. However, to our knowl- 
edge, the full force- velocity relation for ribosomes has not 
been measured so far. But, with the rapid progress in the 
experimental techniques, it should be possible in near fu- 
ture to test the full force-velocity relation predicted by 



I 

our theory. 

We have presented a few quantitative characteristics 
of the fluctuations in the kinetics of ribosomes. We have 
defined a "randomness parameter" r, which is a mea- 
sure of the fluctuations in the dwell times. From the full 
probability density of the dwell times, we have derived 
the expression for r and analyzed some of its interesting 
features. We have also reported the analytical expres- 
sion for the diffusion constant and related it to the mean 
velocity and the randomness parameter. Using the cen- 
tral limit theorem, we have argued that the distribution 
of the run times of the ribosomes from the start codon 
to the stop codon is Gaussian and also pointed out the 
relations between its first two moments and those of the 
dwell time distribution. 

To our knowledge, the run time distribution of ribo- 
somes has not been measured so far. RNA polymerase 
(RNAP) motor runs on a DNA track using the track 
to polymerize the complementary RNA. There are some 
similarities between template-dictated polymerizations 
driven by ribosome and RNAP. The run time distribu- 
tion of RNAP has been measured and found to be Gaus- 
sian [5^. This is consistent with the Gaussian run time 
distribution for ribosomes predicted in this paper which 
follows from very general arguments based on the central 
limit theorem. 

Incorporating inter-ribosome steric interactions in the 
model, we have developed a model for ribosome traffic. 
The model may be regarded as a TASEP for hard rods 
each of which has five distinct "internal states"; tran- 
sitions between these internal states constitute parts of 
the mechano-chemical cycle of a ribosome. Initiation and 
termination of the polymerization of individual proteins 
are captured by imposing open boundary conditions. For 
this model, we have drawn three-dimensional phase di- 
agrams in spaces spanned by parameters which can be 
varied in a controlled manner in laboratory experiments 
in-vitro. In principle, the phase diagram can be obtained 
by analyzing the density profile of the ribosomes in elec- 
tron micrographs of the system for several different con- 
centrations of amino acid subunits, GTP concentration 
etc. 

Appendix I 
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The solution of the equations H])®, for the initial con- 
dition (fTO|) are given by 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 
t 



FIG. 15: (Color online) The probabilities P^{t) {/i = 
1,2, ...,5), corresponding to the initial conditions Pi{0) = 1 
and P2(0) = PsiO) = P4(0) = PsiO) = (see equations f79|l - 
(|83|) ') are plotted for the parameter set uja ~ 23.0, luhi = 24.0, 
k2 = 25.0, LUg = 26.0, LUh2 = 27.0 and cop = 0.5. 



Pl{t) = exp{-UJa't) + WaWp 



exp{—ujit) 



exp{—uj2t) 



exp{~LOat) 



(79) 



exp{—ijJit) exp{~bJ2t) 



(80) 



ea;p(— Wit) 



ea;p(— c<;2i) 



exp{—k2t) 



(CJ2 - ^^l)(fc2 - ^^l) (cJi - LJ2)(fc2 ~ W2) (^1 - fc2)(^^2 - fe) 



(81) 



P4{t) = UJaLLlhlk2 



exp{—LUit) 



exp{~uj2t) 



{UJ2 - Wi)(fc2 - UJi){LLlg - LUi) {lOi - CJ2)(fc2 - W2)(Wg - LU2] 

exp{—k2t) exp{—ujgt) 



(CJI - A;2)(t^2 - fc2)(t^5 - fe) (Wl - Wg)(a;2 - - t^g) 



(82) 



ea;p(— Wit) 



exp{—uj2t) 



{UJ2 - Wl)(fc2 - i^l){0Jg - L0l){'^h2 - ^l) {oJi - LU2){k2 - U}2){^g - 0J2){uJh2 - W2) 

exp{—k2t) ^ exp{~LUgt) 

(Wl - k2){uJ2 - - k2){uJh2 - ^2) (t^l - ^g){^2 ~ LUg){k2 ~ LUg){uJh2 - i^g) 

exp{-UJh2t) 

{UJI - t^ft2)(t^2 - i^h2){k2 - (^h2){(^g - LUh2) 



(83) 



These distributions are plotted in fig[T5] for one set of 
values of the model parameters. These clearly shows that 
the probability P\{t) decreases monotonically from the 
initial value 1 while the states 2, 3, 4 and 5 "rise" and 
"fall" in a sequence. 
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